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We present a lattice QCD calculation of the parameters a and /3 which are necessary in the theoretical estimation 
of the proton lifetime in grand unified theories (GUTs) using chiral lagrangian approach. The simulation is carried 
out using the Wilson quark action at three gauge coupling constants in the quenched approximation. We obtain 
|a(2GeV)|=0.0091(08)(j:}9)GeV^ and |/3(2GeV)|=0.0098(08)(+2o)GeV^ in the continuum limit where the first 
error is statistical and the second one is due to scale setting. 



1. Introduction 

Proton decay is one of the most important pre- 
dictions of grand unified theories (GUTs) and a 
lot of experimental effort has been devoted to de- 
tect it. Up to now, however, no such decay pro- 
cess have been observed, and the lower bound of 
the proton lifetime has been pushed up to 0(10'^^) 
years by the Super-Kamiokande experiment. 

To constrain the parameter space of GUTs by 
the proton lifetime, we need a reliable estima- 
tion of the hadronic matrix elements {PS\0\N), 
where PS and N represents pseudoscaler me- 
son and nucleon respectively and O is a three- 
quark operator. In principle, lattice QCD allows 
a precise determination of the matrix elements 
from the first principles. In previous work, we 
made a GUT-model-independent calculation of 
the matrix elements for all dimension-six three- 
quark operators using the Wilson quark action 
and the plaquette gauge action in the quenched 



approximation Tl . 

There remain two major systematic errors in 
the previous calculation: scaling violation and 
quenching effects. As a first step to reduce the 
systematic errors, this work is devoted to study 
the scaling violation effects. We leave incorpo- 
ration of the dynamical quark effects to future 
work. 

2. Method 

With the aid of chiral perturbation theory [2], 
the proton decay matrix elements {PS\0\N) are 
related to the so-called a and (3 parameters, which 
are defined by 
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where Up denotes the proton spinor. Calcula- 
tion of these parameters is much simpler than 
the direct calculation of the matrix elements 
{PS\0\N): the former is obtained from two-point 



2 



Table 1 

Simulation parameters. The lattice spacing a is 
determined from m„. 



p 


L '-^ xT 


a~'[GeV] 


#conf. (Ref. 3 J 


5.90 


32^x56 


1.934(16) 


300 (800) 


6.10 


40^x70 


2.540(22) 


200 (600) 


6.25 


48^x84 


3.071(34) 


140 (420) 



functions with zero spatial momentum projec- 
tion, while the latter requires three-point func- 
tions with finite spatial momenta. It should be 
noted that our previous study showed that the 
matrix elements obtained from the three-point 
functions are roughly comparable with the tree- 
level predictions of the chiral lagrangian with the 
a and /3 parameters determined on the lattice. 

The lattice operators relevant for the a and f3 
parameters are renormalized as 



(a) 



where 
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and the additional mixing operator is 



O. 



7m J 



(3) 



(4) 



(5) 



The renormalization constants Z, Z^i^, and -^j'^jx 
are calculated perturbatively at one- loop levelpp. 
The continuum operators are defined in NDR 
(naive dimensional regularization) with the MS 
subtraction scheme. 

To obtain the matrix elements, we compose a 
ratio 



R{t) 



(6) 



where ^/Z^ is defined by 
,(^') 



Jp,,(0,0)|p(^')(0)) = 

Zpuf ' . This factor is estimated from the pro- 
ton correlator with local source and local sink. 
It is well known that a precise determination of 
is rather hard because of large statistical 
fluctuations of the local-local correlator. On the 
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Figure 1. Typical light quark mass dependence 
of the a parameter. 

other hand, the ratio of two-point functions can 
be calculated using the smear-local proton corre- 
lator whose statistical fluctuation is fairly small 
compared to the local-local correlator. 

This situation leads us to the following strat- 
egy: (i) For ^/Z^ we resort to the result obtained 
from high statistical calculation of the quenched 
light hadron spectrum performed by the CP- 
PACS collaboration^], (ii) Employing the same 
parameters as the CP-PACS calculation, we make 
a new simulation to estimate the ratio of two- 
point functions including the mixing operators, 
which were not incorporated in the previous CP- 
PACS simulation. After combining the results of 
and the ratio of two-point functions we can 
achieve a few % level of statistical accuracy for 
the matrix elements at each /? and hopping pa- 
rameter. 

3. Simulation details 

Our calculation is carried out in quenched QCD 
with the Wilson quark and the plaquette gauge 
actions. The simulation parameters are given in 
Tablc^ We adopt three values of lattice spacings 
in the range a w 0.1—0.064 fm. The spatial lat- 
tice size is about 3 fm to avoid finite size effects. 
We employ four quark masses corresponding to 
mps/my ~ 0.75—0.5 for each /3. 

For smearing the source, the quark wave func- 
tion in the pion is measured for each hopping 



parameter on gauge configurations fixed to ttie 
Coulomb gauge except for t—0 time slice where 
the wall source is placed^. To estimate y/Z^ 
we first fit the smeared-local proton correlator 
to a single exponential form Z'pexp(—mpt), and 
then the local-local proton correlator is fitted to 
Zpexp(— TOpi) with rrip fixed to the value deter- 
mined by the smeared-local correlator. 

The matrix elements defined on the lattice are 
converted to those in the continuum at /i = 
1/a and the scale is evolved up to /i=2GeV us- 
ing the two-loop renormalization factor in the 
continuumj^. 
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4. Results 

Figure ^ shows the light quark mass depen- 
dence of the a parameter at P — 6.1, which is 
well described by a linear function. Similar quark 
mass dependence is found for the (3 parameter. 
We check that linear and quadratic extrapola- 
tions yield consistent results within error bars in 
the chiral limit at all the lattice spacings. 

In Fig. 121 we plot the parameters a and j3 in 
physical unit as a function of lattice spacing a, 
where we experiment with three choices m^r, nip 
and /tt as physical input to determine the lat- 
tice spacing. In case of the nucleon mass m^v 
as physical input we find little scaling violation, 
which allows us to take the continuum limit of 
the parameters by linear extrapolation. We take 
the extrapolated value as the central one of the 
a and f3 parameters in the continuum limit. On 
the other hand, the parameters show rather large 
scaling violation if we set the lattice spacing by 
the rho meson mass nip or the pion decay con- 
stant fj!-. Although in these cases the values at the 
continuum limit possibly deviate from the result 
obtained with mN as physical input, simple linear 
extrapolation is not reliable to estimate the devi- 
ation. Instead, we estimate this ambiguity from 
the result of quenched light hadron mass spec- 
trum obtained by the CP-PACS collaboration^], 
which shows that the values of and in 
quenched QCD deviate from the experiment by 
about 7% and 10% respectively in the continuum 
limit, once we set the lattice spacing by mp. 



Figure 2. a (open) and /3 (filled) parameters as 
a function of lattice spacing. 

Taking account of this ambiguity, we obtain 

|a(2GeV)| = 0.0091(08) (1^°) GeV^ (7) 

|/3(2GeV)| = 0.0098(08) GeV^ (8) 

with a/ (3 < in the contimmm limit, where the 
first error is statistical and the second one is due 
to scale setting. We note that the relative sign of 
a and /3 could be important, whereas the overall 
one is irrelevant. 
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